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1. Introduction . In this note we consider the analysis of 
a plate which consists of parallel load-carrying strips of high 
modulus separated by buffer strips of relatively low modulus [1]. 
Both strips are boron-epoxy composite. Through an appropriate 
orientation of the fibers, a relatively low modulus and high 
toughness is achieved in the buffer strips. This in turn improves 
the crack arrest character! sties of the structure as a whole. 

Even though the primary (load-carrying) laminates and the buffer 
strips are highly anisotropic, in this note for reasons of analyt- 
ical expediency, it will be assumed that both materials are iso- 
tropic and linearly elastic (Figure 1). A symmetri cal ly located 
through crack is assumed to form and propagate in the primary 
laminate, and eventually to enter into the buffer strips. The 
objective of this note is to calculate an upper and a lower bound 
for the stress intensity factor when the crack tips are in the 
primary laminate, in the buffer strips, and at the interfaces. 

2. Bounds for the stress intensity factor . The actual 
problem is described in Figure 1. The composite plate is assumed 
to contain a symmetrically located through crack of length 2a, 

*This work was supported by NASA-Langley under the Grant 
NGR 39-007-011. 
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where 0<a< h 2 , h 2 = h ] + 2b 2 , 2h-j and 2b 2 being the widths of 
the primary laminate (material 1) and the buffer strip (material 
2), respec ti vely . The crack is perpendi cul ar to the bimaterial 
interfaces. The plate is subjected to uniaxial loading parallel 
to the strips outside the perturbation zone of the crack. The 
elastic constants of the materials 1 and 2 are E-j , (or 
]j-j - E i / 2 ( 1 +v-| ) , K-| = (3~v-j)/(l^v^)j and E 2 , ^2 P * ^2^* 

respectively. The elastic constants of the composite plate may 
be evaluated from 

E c = ( h-j E-j + b ^ E 2 ) / (^i h 2 ) * 

v c = ^ h l v l + b 2 v 2 ^ ^ ( b i + b 2 ) . (1) 

In the plate without the crack the stresses o-j and a 2 acting on 
the strips in a plane perpendicular to the interfaces are related 
by 

a 1 /a 2 = E ] /E 2 . (2) 

The solution of the problem described in Figure 1 is not 
available and, at present, appears to be intractable. However, 
through a judicious choice of the elastic constants and the 
dimensions in the problem described in Figure 2 (the solution of 

which is available [2]) one could find estimates for an upper and 

a lower bound for the corresponding stress intensity factors. 

One simple procedure to obtain these estimates is the following: 

(a) An upper bound for the stress intensity factor k when 
the crack is in the primary laminate: 0< a < 

This estimate may easily be obtained by replacing the 
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part of the composite plate |x[ > h 2 (see Figure 1) by homogeneous 
half planes having the elastic constants of the buffer strip, 

^2* v 2* Since the stiffness of the resulting half planes |x| > h-j 
is less than that of the half planes in the actual structure, the 
stress intensity factor calculated from Figure 2 with 

p ( x ) “”0^, a = a £ h-j , h ** h -j , E-j ~ E ^ , 

v* = v, , l* z = E 2 , v* = v 2 , (3) 

would be greater than the k value corresponding to Figure 1. 


(b) A lower bound for k when the crack is in the primary 
laminate: 0 < a < h-j . 

This estimate may be obtained by replacing the non- 
homogeneous composite plate for |x|>h-j by the equivalent homo- 
geneous composite half planes where 


p(x) = -a 1 


= a < h| , 


= h 1 , 


E ! " E l’ 


rk 

V-j = V-| 


F* = F 

t c 


★ 

v 2 = v c 


(4) 


Here the reason for the resulting k being lower than the actual 
k is that the stiffness E c used for calculations in Figure 2 is 
greater than E^ which is the stiffness of the neighboring media 
in the actual problem (Figure 1). 


{ c ) A lower bound for k when the c rack is in the buffer 
strip: < a ;< h 2 . 

This estimate may be obtained by replacing the nonhomo- 
geneous composite plate by a homogeneous strip ( E 2 * v 2 ) (which 
has a stiffness smaller than the actual) for |x|<h 2 and by 
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homogeneous half planes (E-j , v-| } (which have a stiffness greater 
than the actual) for |x| > h 2 and using the pressure p=-o 2 the 

external load. Thus in Figure 2 the parameters would be 

p ( x ) = -a 2 , a * = a » h* = h 2 > E* = E 2 * 

v* = v 2 , E* - E 1# v* = v r (5) 


(d) A closer lower bound for k when the crack is in the 
buffer strip: h-j < a < h 2 . 


This estimate is obtained by using the same plate as 
in (c) with the actual (and higher) crack surface pressure, i.e., 

a = a , h = h 2 , E ~ E 2 , v -j = v 2 , E 2 = E -j , 





for ] x | < h-jV 
for h-j < | x | < a . 


( 6 ) 


Theoretically, the k value given by this estimate is still a 
lower bound, but is believed to be quite close to the actual 
stress intensity factor. 


In the problems described above the solution is obtained 
through a superposition where the stress intensity factor is 
calculated from the singular solution in which the crack surface 
pressure is the only external load (Figure 2). In the solution 
given in this note, as long as the crack tip is away from the 
interfaces, the following asymptotic standard expression is valid 
for the cleavage stress a around the crack tip: 

J V 

o ( x , 0) * • - — , (a/h., 1 = 1,2). (7) 

/2 ( x - a ) 
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However, if the crack tip terminates at an interface, the asymp- 
totic stresses around the crack tip are given by [3,4] 

a _•• ( r , e ) = — — f i n - ( e ) + 0(r a2 ) , (i,j = r,0; 

Re(a 2 )>0, 0 < a < 1 ) , (8) 

where the crack is assumed to lie along 0 = tt (i.e., the negative 
x axis) [3,4]. Foil owi ng Figure 3, if p 2 > p-j then 0 < a < 0 . 5 and 
if < ^1 then 0.5<a<l. The functions f-jj(e) giving the angu- 
lar variation of the stresses around the crack tip for a symmet- 
rically loaded crack (i.e., first mode) and for various values 
of p 2 /u-j are s ^ own figures 3-5 in the following normalized 
form: 

f ij (e) = G 1j (e)/G ee (0) , (i,j = r,6). (9) 

3. The Numerical Results .' The numerical results given in 
this section are obtained for the boron-epoxy composite plate 
which was considered in [1]. The material constants of the 
(equivalent) isotropic materials are: 

Primary laminates: E-j=11.39xlo^psi, v-j = 0.33 , 

Buffer strips: E 2 = 2 . 42x1 O^psi , v 2 =0.85. 

Composite plate (h-^/b-j =2): E c = 8.4xlo^psi, v £ = 0.503, 

(i.e., u-j/]i 2 =6.55, k-j = 2 , k 2 - 1 . 1 6 , p c /u-j = 0.657, 

k s 1,66), 

v 

The results corresponding to the cases (a).-(d) described in 
Section 2 are given in Tables 1-4, respectively, and are summa- 
rized in Figure 6. Note that Table 1 is valid for any value of 
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Figure 3. Angular variation of o 00 around a crack tip 
touching the interface. 



Figure 4. Angular variation of o , around a crack tip 
touching the interface. 



Figure 5. Angular variation of a around a crack tip 
touching the interface. r 


5.0 


—, 5.0 



Figure 6. Bounds for the stress intensity factor in the 
composite plate. 
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the width ratio h -| / b 2 and for 0<a<_h*|. Even though the results 
given in Table 3 are obtained only for h 2 / b -j = 2, it is clear that 
since the pressure on the crack surface is uniform, the results 
for any other value of l^/h^ can be obtained from those given for 
h 2 / h i = 2 by using appropriate a/h^ ratios. Thus Figure 6 shows 
the curves for four different values of b^/h-j ratio, namely 2, 
1.75, 1.5, and 1.25. 


Table 1. Upper bound for the stress intensity factor: 0< a <_ h -j 9 

(crack in primary laminate), k-j = 2, k2 s 1.1*>» ^2^^! = 0-1527 , 
uniform crack surface pressure, -a-j . 


a/h-j 

0.14 

0.225 

0.325 

0.6 

0.75 

0.9 

1.0 

k / a -j /E 

1.007 

1 .018 

1.038 

1.157 

1 .296 

1.618 

(0.550) 
a = 0.66980 


Table 2. Lower bound for k: 0 < a ;< h“| , k-j = 2 , k2 = k c =1.66, 

y c /p'j = 0. 657 , h -j / b 2 - 2, (crack in primary laminate), uniform 
crack surface pressure, -a-j . 


a/h 1 

0.6 

0.75 

0.9 

1.0 

k / a 1 /a 

1 .031 

1.055 

■ 

- — - ■ — 1 

1.099 

(0.920) 
a = 0.52620 


Table 3. Lower bound for k: hg>_a> h-|» (crack in buffer strip), 

k-| = 2 , K 2 = 1 ‘ 16 > uniform crack surface pressure , - o 2 - 


a/h 2 

0.5 

0.75 

0.875 

1.0 

k/Og/a 

0.899 

- 

0.772 
' 

0.675 

(1.160) 
a = 0.42724 
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Table 4. A (closer) lower bound for k: h 2 ^a>h] (crack in 

buffer strip), kj = 2, *2 = 1.16, ij-|/p2=6.55, crack surface 
pressure p(x) = for 0 <_ | x | £ h-j , p(x ) = -og f° r h-j < l x l <a > 
P2/ P] " f^/E-j * 


h l . 1 
h 2 ' 1 

a/h 2 

0.5 

0.625 

0.75 

0.875 

1.0 

k/c 2 /a 

4.231 

2.396 

1 .820 

1.302 

(2.221 ) 
a = 0.42724 

H 1 1 
h 2 T775 

a/h 2 

1/1.75 

0.68 

0.92 

1 .0 



k/a 2 /a 

4.120 

2.580 

1.397 

(2.200) 
a = 0.42724 

h l 1 

h 2 - 77T 

a/h 2 

2/3 

0.8 

0.94 

1.0 

k/o 2 /a 

3.85 

2.257 

1 ,536 

(2.900) 
a = 0.42724 

h l 1 

a/h 2 

0.8 

0.88 

0.96 

1.0 

h 2 “ 1775 

k/a 2 /a 

3.43 

2.286 

1 .666 

(3.650) 
a = 0.42724 


Table 4 shows the results corresponding to the assumption ( d ) 
of the previous section for four different values of h 2 /h-|. It 
should again be noted that these results give a (closer) lower 
bound for k when the crack is in the buffer strip. Since in this 
case changing the material constants to obtain an upper bound is 
not realistic, an approximation to an upper bound may be obtained 
by changing the crack surface tractions. This may be done by 
using the highest possible uniform crack surface pressure, p = ra-j 
instead of the lowest uniform pressure p = -c? 2 1n t ^ ie calculations. 
The result would be simply the multiplication of the values found 
in Table 3 and Figure 6 (c) by the factor o-j/c^ = = 4.7066 
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(see, for example, the dashed curve (c 1 ) in Figure 6 which corre- 
sponds to i^/h-] =2). These results are repeated in Table 5. 


Table 5. An (approximate) upper bound for k: (crack 

in buffer strip), k-| = 2 , < 2 = 1 . 16 , p (x ) = -aj , 

(0 <J x | < a) . 


a/h 2 

0.5 

0.75 

0.875 

1.0 

k / a 2 ^' 

4.231 

3.634 

3.177 

(5.460) 
a = 0.42724 


The tables 1-5 also give the results for the crack terminat- 
ing at a bi -material interface (i.e. , a = h^ or a = h^). The 
corresponding k values shown in the tables in parentheses and the 
exponents a should be understood (and used) in the sense as 
defined by eq. ( 8 ). Note that the k values given in the tables 
and in Figure 6 tend to infinity or zero as the crack tip 
approaches the interface. The reason for this, of course, is 
found in the definition of k as given by eq . (7) and in the 
abrupt change in the power of the singularity a when the crack 
tip touches the interface. A possible technique of applying 
these results regarding the irregular singularity a^ 0.5 in 
fracture studies may be found in [5]. 
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